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Abstract 

The time-dependent response of polystyrene and poly(methyl methacrylate) is 
studied in isothermal long-term shear creep tests at small strains and various tem- 
peratures in the vicinity of the glass transition point. A micromechanical model 
is derived to describe the experimental results. Constitutive equations are devel- 
oped under the assumption that the behavior of amorphous polymers is governed by 
two micro-mechanisms: rearrangement of cooperatively relaxing regions (CRR) re- 
flects the viscoelastic response, whereas displacement of CRRs with respect to each 
other is responsible for the anelastic response. It is demonstrated that some crit- 
ical temperature exists slightly above the glass transition temperature, where the 
dependences of adjustable parameters on temperature are dramatically changed. 
The critical temperature is associated with transition from dynamic heterogeneity 
in amorphous polymers to static inhomogeneity. 
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1 Introduction 



The mechanical behavior of amorphous polymers in the vicinity of the glass transition 
point, Tg, has attracted essential attention in the past decade, see, e.g., monographs 
The researches focus conventionally on three main issues: 

1. the kinetics of slow out-of-equilibrium evolution of internal structure for polymers 
equilibrated at some temperature above Tg and quenched to a temperature T in the 
sub-Tg region (physical aging of polymeric glasses H); 

2. the existence of some critical temperatures, Tc, slightly above or below the glass 
transition point, associated with the crossover point predicted by the mode-coupling 
theory 0, as well as with splitting points for a and slow /? relaxations [0, |l^ 



and for a and a relaxations |I3| ; 

temperature-induced changes in the mechanical response of amorphous polymers 
from the "glassy-like behavior" in the sub-Tg region to the "rubbery-like behavior" 



in the post-Tg domain observed in standard creep and relaxation tests [[1^] . 

The present paper is concerned with the latter two questions. 

A conventional standpoint is that the entire information about the time-dependent 
behavior of solid polymers at small strains can be obtained in isothermal short-term tests 
(with the duration of about 10^ s) and the mechanical response in creep and relaxation 
tests can be adequately predicted in the framework of linear viscoelasticity. To demon- 
strate that this is not the case, we perform long-term creep tests (with the duration of 10^ 
s) on two widely used amorphous polymers, polystyrene (PS) and poly(methyl methacry- 
late) (PMMA), at various temperatures below and above Tg. Experimental data reveal 
rather poor superposition of creep compliance curves plotted in double logarithmic co- 
ordinates by their shifts along the time-axis: superposition of the middle parts of the 
long-term compliance curves results in deviations of about 10 to 15 % (in the logarithmic 
scale) in the regions of small and large times. This observation implies the conclusion 
that viscoelasticity is not the only mechanism that governs the time-dependent behavior 
of polymers at the macro-level. To describe another possible mechanism, the following 
scenario is proposed for deformation of amorphous polymers: 

1. At some temperature Tf > Tg a homogeneous state of an amorphous polymer be- 
comes thermodynamically unstable, and dynamic micro-heterogeneities arise in the 
bulk material. These inhomogeneities (Fisher's clusters |]15| or frustration-limited 



domains |jT6|) are observed as temporary micro- regions with higher mass densities 
compared to the density of the surrounding domains. The presense of the Fisher 
clusters is confirmed experimentally in light scattering tests which reveal long-range 
density fluctuations (with the lengthscale of 200 nm) in molecular and polymeric 
liquids slightly above the glass transition temperature ITT 



2. Below Tp, an amorphous polymer is treated as an ensemble of high-density clusters 
bridged by long chains (which create "less cohesive spaces" between "more cohesive 
domains"). This picture is confirmed by low- frequency Raman spectroscopy for 
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poly(methyl methacrylate) []TB| and polycarbonate |]T^ and by X-ray diffraction and 
electron microscopy for poly(ethylene tereplithalate) [|20[. Observations evidence 
that the less cohesive domains create a penetrating network, whereas more cohesive 
regions may be thought of as isolated islands with the characteristic size of 4 to 7 
nm (below the glass transition point). The latter result is in quantitative agreement 
with the characteristic lengths of cooperative regions determined for ortho-terphenyl 
salol confined to nanopores p2[ and mixture of polychlorinated biphenyls 



by means of dynamic light scattering, photon correlation spectroscopy and dielectric 
relaxation. The islands are separated by the less cohesive space with the average 
distance between high density clusters of about 1 nm [jl8| . 



The domains with high density are identified with cooperatively rearranged regions 
(CRR) predicted by the concept of cooperative relaxation Any CRR consists 



of dozens of strands linked by entanglements and van der Waals forces |25] and 
may be treated as a "ball-like structure" with a paracrystalline-type order . Re- 
arrangements of CRRs occur at random times as they are thermally agitated. A 
rearrangement event is treated as large-angle rotation of segments of strands in a 



CRR [^. Large-angle reorientation for molecular clusters was recently confirmed 
by multidimensional deuterium nuclear magnetic resonanse spectroscopy for super- 
cooled glycerol and toluene ||28|. Rearrangement of CRRs at the micro- level is 



associated with the viscoelastic response of a specimen at the macro-level. 

4. The less cohesive space between CRRs has a complicated fractal structure. It is 
modeled as an ensemble of mutually independent less cohesive domains (LCD) with 
the characteristic size of about 100 nm (the length scale of long-range density fluc- 
tuations observed in light scattering experiments) . Two main functions are ascribed 
to LCDs: (i) to transmit macro-strains in a sample to individual CRRs, and (ii) 
to ensure sufficient freedom for mutual displacement of CRRs with respect to each 
other (within a LCD to which they belong). Micro-deformations driven by dis- 
placements of CRRs are associated with the anelastic response of a specimen at 
the macro-level. Transformation of the internal structure of less cohesive space 
driven by mutual displacements of CRRs was confirmed by low-frequency Raman 
spectroscopy for poly(methyl methacrylate) at finite shear deformations |^ . 



The meaning of the term "anelastic deformation" used in this study slightly differs from 
that conventionally employed in the analysis of the viscoelastic and viscoplastic behavior 
of amorphous polymers far below their glass transition point, see, e.g., |2^, According 
to the concept of shear micro-domains |Q , the anelastic deformation is treated as a defor- 
mation which recovers during the experimental time scale at temperatures below Tg — 20 
K, whereas the plastic deformation is thought of as that recoverable in the close vicinity of 
the glass transition temperature only [0. To establish some correspondence between our 



definition and the traditional one, we treat the less cohesive space as a rubber-like material 
and assume that some chains in the less cohevise domains are disentangled when rela- 
tively large stresses are applied to an amorphous polymer. The disentanglement results 
in a substantial increase in the average size of a LCD driven by coalescence of neigh- 
boring LCDs (which receive sufficient freedom to be aggregated). This implies that the 
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mechanical responses of a sample at loading and subsequent unloading do not coincide, 
which is evidenced as the onset of residual strains. These strains disappear only when 
a sample is heated to the glass transition region and annealed for a time sufficient for 
healing (reformation) of entanglements (which is tantamount to the return of the internal 
structure of a polymer to its virgin state). 

It is worth noting that a similar process (disentanglement of chains in the less cohesive 
space and aggregation of LCDs into larger clusters) can occur with an increase in temper- 
ature above the glass transition point. In this case, however, thermal fluctuations may be 
thought of as a driving force for the disentanglement process instead of mechanical stresses 
(because the intensity of fluctuations grows with temperature, T). This implies some sim- 
ilarity between viscoplastic deformations of amorphous polymers below Tg and their flows 
at temperatures far above Tg. This similarity was recently conflrmed by comparison of 
the activation enthalpies for plastic deformation of glassy polycarbonate, poly(methyl 
methacrylate), polystyrene, poly (vinyl chloride) and styrene-acrylonitrile copolymer with 
those for the non- Newtonian flow of their melts BTI, 15^ . 



The objective of this study is three-fold: (i) to analyze the response of PMMA and 
PS in isothermal long-term creep tests with small strains at various temperatures, (ii) 
to derive constitutive equations for the behavior of amorphous polymers based on the 
two micro-mechanisms for the time-dependent response (rearrangement of CRRs and 
displacement of CRRs with respect to each other) and (iii) to demonstrate the presence 
of some critical temperature, Tc, slightly above the glass transition point, where the 
dependences of adjustable parameters on temperature dramatically change. 

The exposition is organized as follows. Section 2 is concerned with the description 
of observations. Constitutive equations for the time-dependent response of amorphous 
polymers are developed in Section 3. In Section 4, material constants in the stress-strain 
relations are found by fltting experimental data. Section 5 deals with a discussion of our 
findings. Some concluding remarks are formulated in Section 6. 



2 Experimental 

The materials studied are commercially available polystyrene (the weight-average molecu- 
lar weight Myj = 352 kg/mol and the polydispersity M^/Mn = 2.2) produced by Hoechst 
AG and poly(methyl methacrylate) [My^ = 98 kg/mol, Myj/Mn = 1.6) produced by Rohm. 
To determine the specific volume of polymers, v, as a function of temperature, T, volume 



dilatometric measurements were carried out following the procedure exposed in [Q. A 
specimen with the mass 3 g was held at the initial temperature Tq = 190 °C for 30 min 
to erase the effects of thermo-mechanical prehistory, and, afterwards, was cooled with the 
rate 3 K/h to the final temperature = —10 °C. The error in dilatometric measurements 
does not exceed 10~^ cm^/g. 

The specific volumes of PS and PMMA are plotted versus temperature in Figures 1 
and 2. Experimental data are approximated by the linear dependence 

V = Vo{l + aoT), (1) 

where the adjustable parameters vq and ao are found by the least-squares technique 
separately in the regions of low (curve 1) and high (curve 2) temperatures. The glass 
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transition temperature, Tg, is determined as a point of intersection for straight lines 1 
and 2. The values of Tg and the coefficient 

dv 

a = — = aovo 
ai 

are listed in Table 1. The results are in quantitative agreement with data provided by 
other sources (see, e.g., Tg = 100 for PS and Tg = 105 °C for PMMA in [I and 91.2 
for PS and 98.7 °C for PMMA in p3[). The difference between our experiments and the 



33| , consists in the following: 



1. in the present study polystyrene specimens were preliminary melted in vacuum to 
remove additives (paraffinic wax), whereas Greiner and Schwarzl used them as 
received; 

2. different commercial grades of poly (methyl methacrylate) were tested. 

Torsional creep was measured at the stress levels that ensure the linear response of 
polymers by using a high-accuracy creep tester which was previously described in detail 
in . The samples have the shape of circular cylinders with height 105 mm and radius 6 



mm. The specimens were heated to the temperatures To = 115 °C (PS) and To = 130 °C 
(PMMA) and held at this temperature for at least 30 min to erase the effects of thermo- 
mechanical prehistory. Any specimen was quenched to the test temperature, T, in the 
apparatus. After the preconditioning time of at least 1 h, the specimen was loaded with 
a constant stress, and the torsional creep compliance, J, was measured as a function of 
time t. The response of PS at temperatures above 115 °C was measured after heating of 
samples in the apparatus to the test temperature, T, and preconditioning for 1 h. The 
response of PMMA below the glass transition temperature was measured on specimens 
which were partially equilibrated (for this purpose, the preconditioning time was increased 
from 1 to 155 h depending on the test temperature). 

Figures 3 and 4 show the creep compliances of PS and PMMA at various temperatures 
in the range from 95 to 140 °C for PS and from 92.5 to 130 °C for PMMA, which, according 
to Table 1, correspond to the region of rubbery response for PS and to the regions of glassy 
and rubbery behavior for PMMA. 



3 Constitutive equations 

An amorphous polymer is treated as an ensemble of cooperative rearranging regions 
bridged by less cohesive space. Denote by Xq = Xq(T) the number of CRRs (per unit 
mass) at temperature T and by m = m{T) the average mass of a relaxing region. The 
total mass of CRRs occupying unit mass of the polymer reads Mq = mXo, and the total 
mass of LCDs (per unit mass) is given by 

M = 1 - Mo = 1 - mXo. (2) 

The macro-strain in a specimen, e, is assumed to coincide with the sum of the micro-strain 
in CRRs and LCDs, e, (for simplicity, we confine ourselves to uniaxial deformation and 
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assume that micro-strains are the same for all CRRs and LCDs) and the micro-strain, Ca, 
induced by displacement of CRRs as rigid bodies with respect to each other, 



e{t)^e{t) + e,{t). (3) 

The subscript index "a" indicates that mutual displacements of CRRs at the micro-level 
are associated with the anelastic response of a polymer at the macro-level. 

At small strains, a LCD is modeled as a linear elastic medium with the strain energy 
density (per unit mass) 

where the rigidity of CLD per unit mass of a polymer, //lcd, is expressed in terms of its 
rigidity per unit volume, clcd, by the formula 

M 

A^LCD = Clcd — , 
P 

where p is mass density of LCD (we suppose that the difference between the densities of 
more cohesive regions (CRR) and less cohesive space (LCD) is small, which implies that 
p may be associated with the mass density of a polymer). 

To calculate the mechanical energy of CRRs, we introduce the function X{t, r) which 
equals the number of CRRs (per unit mass of the polymer) at time t that have last 
been rearranged before instant r. The function X entirely determines the rearrangement 
process: 

• X{0, 0) coincides with the concentration of CRRs in a stress- free medium, 

X(0,0) = Xo, (4) 

• the expression 



dX ^ , 



dr 

t=T 



equals the number of CRRs (per unit mass) rearranged within the interval [r, T+dr], 
the quantity 

-^{t,r)dT 

equals the number of these CRRs that have not been rearranged during the interval 
the amount 

dX , , , 
- — (t,0)dt 

equals the number of CRRs (per unit mass) that have been rearranged for the first 
time during the interval [t, t + dt], 
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• the expression 

determines the number of CRRs (per unit mass) that were rearranged within the 
interval [r, r + dr] for the last time before their rearrangement during the interval 
[t,t + dt]. 

The viscoelastic response of a polymer is modeled as a sequence of (driven by thermal 
fluctuations) random hops of rearranging regions in their potential wells PH|. According 



to the trapping concept any CRR is modeled in the phase space as a material 

point trapped at the bottom level of its potential well on the energy landscape. At 
random times, the point hops to higher energy levels as the CRR is thermally agitated, 



but it cannot leave its trap. Unlike previous studies, see, e.g., |^|, all potential wells are 
assumed to have the same depth. This allows the number of adjustable parameters in the 
model to be substantially reduced. The account for the distribution of energies for cages 
does not lead to an essential improvement in the quality of matching experimental data 
and does not result in qualitative changes in our conclusions. 

Referring to the transition-state theory we postulate that some liquid-like (refer- 
ence) energy level exists on the energy landscape, where CRRs change their configurations. 
When a CRR reaches the liquid-like state in a hop, stresses totally relax in it. If a relax- 
ing region hops below the reference level, it returns to its position at the bottom of the 
potential well without changes. 

Let us consider a CRR at time t that has last been rearranged at instant r G [0, t) and 
denote by eo(t, t) the strain from its stress-free configuration at time r to the deformed 
configuration at time t. Because a CRR totally relaxes reaching the liquid-like state, its 
stress-free configuration coincides with the deformed configuration of the bulk material 
at the instant of rearrangement. This implies that the strain eo(t, r) is given by 

eo(t,r) =e(t) -e(r). (5) 

At small strains, a CRR is thought of as a linear elastic medium with the mechanical 
energy 

ucRR{t,T) = -ceo(t,r), 

where c is the rigidity per CRR. Multiplying the mechanical energy, mcrr(^, 0), of CRRs 
that have not been rearranged until time t by their number, X{t,0), and using Eq. (5), 
we obtain the mechanical energy of non-rearranged CRRs (per unit mass). 

Multiplying the mechanical energy ucRR{t,T), of CRRs that have last been rearranged 
during the interval [r, r + dr] by their number dX/dT{t,T)dT, we find the mechanical 
energy of these CRRs, 

cdX 



dr. 
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Summing these expressions for various r e [0,t), we arrive at the formula for the strain 
energy density of relaxing regions, 



UcRRit) = - 



Xit,0)e\t) + ^(t,r){eit)-eir)) dr 



Neglecting the energy of mutual interaction between CRRs and between CRRs and less 
cohesive space, we determine the mechanical energy of a polymer (per unit mass), U, as 
the sum of the strain energy densities for CRRs, Ucrr, and LCDs, C/lcd, 



U(t) ^ - 



clcdM 
P 



+ cX{t, 0) 



dr. 



(6) 



It follows from Eqs. (3) and (6) that the function U is given by 



Uit) 



ClcdM 
P 

c ftdX 



+ cX{t,0) 



eit) - e,{t) 



dr. 



(7) 



The stress, is expressed in terms of the macro-strain, e(i), by means of the conven- 
tional relation 

Substitution of Eq. (7) into this formula implies that 

a{t) = 

- (e,(i) - ea(T)) 



[ci^cdM + cpX{t,0)][e{t) - e^{t) 
rt dX 



dr. 



(8) 



Bearing in mind that 



X{t,0) + l^^{t,r)dr^X{t,t) 



and using Eq. (4) and the conservation law for the number of CRRs per unit mass, we 
transform Eq. (8) as follows: 

/■* dX 

<y{t) = (clcdM + cpXo)[e(i) - e^{t)] - pc —{t,T)[e{T) - e^{T)]dr. 
It follows from this formula and Eq. (2) that 



where 



a{t) = G[[e{t) - eM ' ^ ^M[e{T) - e,(r)]rfr}, 

, cpXo 



G = clcd(1 - mXo) + cpX^, 



(9) 



(10) 



clcd(1 - rnXo) + cpXg ' 

The constitutive equation (9) is determined by two functions, X{t,T) and ea(t), which 
will be determined later. 
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3.1 The kinetics of rearrangement 

Denote by 7 the attempt rate (the average number of hops in a trap per unit time) and 
by Q the probabihty of reaching the hquid-hke energy level in a hop. Multiplying the 
attempt rate, 7, by the probability, Q, we find the rate of rearrangement F, 

r = 7g. 

Equating V to the relative rates of changes in the concentrations of CRRs rearranged at 
various instants, we obtain 



dX 



(t,0) 



ax 



1-1 f)'i 



X(t,0) dt 

The solution of Eq. (11) with the initial condition (4) reads 

dX 



dtdr 



where 



X{t,0) = Xoexp(-n), ^(^'^) = F{r)eM-^{t - r)]. 



^/ X dX^ , 



(11) 



(12) 



T=t 



is the number of CRRs (per unit mass) rearranged per unit time at instant t. To determine 
the value of F{t), we sum the number of CRRs rearranged for the first time at instant t, 

and the numbers of CRRs rearranged (for the last time before t) at some instant r e (0, t) 
and reaching the liquid-like level at time t 

This results in the equahty 



dX , , d'^X , 



dt Jo dtdr 

Substitution of expressions (12) into this formula implies the linear integral equation for 
the function F{t), 

rt 



Fit) = r 



Xocxp(-n) + F{t) exp(-r(t - T))dT 



(13) 



The solution of Eq. (13) is given by 

F(t) = rxo. 

Substitution of this expression into Eq. (12) results in 

dX 

— (t,r) = rXoexp[-r(t-T)]. 
Combining Eqs. (9) and (14), we find that 

a{t) = G|[e(t) - e,{t)] - bV j\xp[-r{t - r))[e(r) - ea(r)]dr|. 



(14) 



(15) 
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3.2 The kinetics of anelastic flow 



To describe the anelastic flow (which is associated with the relative displacement of CRRs 
at the micro- level), we adopt the first order kinetic relation 

de 

V^ = e- e., (16) 

which may be treated as a simplifies version of the flow law for a suspension of rigid 
particles in a fluid, where the left-hand side determines the drag force, whereas the right- 
hand side determines the force acting on the particles from the fluid flow. The parameter 
rj stands for the drag coefficient, which, in general, depends on the deformation history. 
The simplest version of the relation between t] and e (which should be independent of the 
sign of e) reads 

r/ = r/o(l + «:e2), (17) 



where t]q and n are material constants. A similar equality was proposed in where 
the physical meaning of the quantities tjq and k, was discussed in terms of the molecular 
structure of polymers. Substitution of Eq. (17) into Eq. (16) implies that 

where 

^0 

Given a loading program, e = e(t), Eqs. (15) and (18) entirely determine the stress cr(t). 

3.3 Transformation of the governing equations 

Introducing the notation 

S(t) = Tj\xp[-T{t - r))[e(r) - e,{r)]dT, (19) 

we re-write Eq. (15) in the form 

a{t) = G [e{t) - ea(t) - 6S(t)] . (20) 
It follows from Eq. (19) that the function S(t) obeys the differential equation 

^(t) = r[6(t)-e,(t)-S(t)], S(0)=0. (21) 

Nonlinear ordinary differential equations (18) and (21) together with the linear algebraic 
equation (20) describe the viscoelastic and anelastic response of an amorphous polymer. 
In the sequel, we focus on the study of creep tests with 

a{t) = (To (t > 0), 
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where the stress cxo ensures small deformation of a specimen. Setting 



J 



we present Eq. (20) as 



J. 



0-0 



0-0 



Jo 



G" 



(22) 



J{t)^ Jo + Ut) + bS{t). (23) 
Substitution of expressions (22) into Eqs. (18) and (21) results in the differential equations 



l + CJ2(i) ' 

j{t) - ut) - s{t) 



S{0) = 0, 



(24) 
(25) 



where C = nag. Excluding the function S from Eqs. (23) and (25), we arrive at the 
differential equation for the function J, 



dJ. 



^-{t) = ^{t) + r[jo - B(j{t) - J,{t 
where, according to Eq. (10), 



J(0) = Jo, 



1 



1 + r' 



cpXo 



clcd(1 - mXo) 



(26) 



(27) 



Equations (24) and (26) will be used in Section 4 for the analysis of experimental data 
in creep tests. These formulas contain 5 adjustable parameters. A, B, C, Jq and F, 
to be found by fitting observations. The material constants have transparent physical 
meaning: Jq is the creep compliance at the instant of application of external load, F is 
the characteristic rate of relaxation, A is the rate of anelastic flow at the initial interval 
of loading (when an increase in r] driven by the growth of strains is negligible), B is 
the ratio of the stiffness of less cohesive space to the total stiffness of a polymer, and C 
characterizes changes in the rate of anelastic fiow induced by an increase in strains. 



4 Validation of the model 

To find adjustable parameters in the model, we match experimental data depicted in 

Figures 3 and 4. For a creep curve at any temperature, we find the material constants Jq, 
F and B which minimize the discrepancies between observations and results of numerical 
simulation at the first several decades of time (when logt changes from —1 to 3 for PS 
and from to 3 for PMMA). The parameters Jq, F and B are determined by using the 
steepest-descent procedure (to simplify the analysis, A is set to be zero at this stage of the 
analysis). Afterwards, we fix the values of these constants found by fitting experimental 
data in short-term creep tests, and determine the remaining two parameters, A and C, 
which provide the best approximation of observations when logt changes from 3 to 6. 
The opportunity to divide the entire interval of measurements into two sub-intervals 
(short- and long-term data) is based on the observation that the effect of anelastic flow 
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(which is described by the parameters A and C) is neghgible during the short-term tests, 
whereas the influence of the viscoelastic response (which is described by the constants Jo, 
r and B) becomes rather small when time, t, essentially exceeds the characteristic time of 
relaxation, F"^. Figures 3 and 4 demonstrate fair agreement between results of numerical 
simulation and experimental data for the two amorphous polymers. 

To evaluate the effect of temperature on the material constants, we plot the parameters 
F, A, B, C and G = Jq^ versus T. Because the viscoelastic behavior is associated at the 
micro-level with rearrangement of CRRs driven by thermal fluctuations, the dependence 
F = F(T) is assumed to be described by the Arrhenius law 

F = F.exp(-^ 

where F* the rate of thermal fluctuations at high temperatures, R is the universal gas 
constant and AE-p is the energy of activation. It follows from this equality that 

logF = Fo-^, (28) 



where 



AEr 

log F„ Fi = — ^ log e, log = logio 



Figures 5 and 6 evidence that Eq. (28) describes experimental data fairly well when the 
adjustable parameters Fj are determined independently in three different regions: 

1. the low-temperature region below the glass transition point; 

2. the medium-temperature region between the glass transition temperature, Tg, and 
some critical temperature, T^, 

3. the high-temperature region where temperature, T, exceeds T^,. 

The values of Tg and found by matching observations are collected in Table 2. 

To compare our results with those available in the literature, we calculate the fragility 



index m, see ^9[, which, in our notation, is defined as 

rflnF 



m 



T=T„ 



d{TjT) 

It follows from this equality and Eq. (28) that 

m = ^lnlO. (29) 

We calculate rh for the two polymers by using Eq. (29) and the data depicted in Figures 5 
and 6 and find that m = 121.2 for PS and m = 98.9 for PMMA. These quantities provide 
excellent agreement with the data collected in |^^: rh = 116 for PS and rh = 103 for 
PMMA. 
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The anelastic flow of CRRs is assumed to be tliermally activated, wliicli implies that 
the dependence A — A{T) may also be approximated by the Arrhenius law with some 
characteristic rate and activation energy AEa- By analogy with Eq. (28), we write 



1 A A 

logA^Ao- — 



(30) 



with 



Ao = log A 



A, 



AEa 



loge. 



R 



Experimental data for A{T) arc depicted in Figures 7 and 8 which show that Eq. (30) 
provides an acceptable approximation of observations. 

The adjustable parameters B, C and Go are plotted versus temperature T in Figures 9 
to 14. The effect of temperature is fairly well described by the phenomenological relations 



where AT — T — Tg and the constants Bi, Ci and are determined by the least- 
squares technique. The quantity AT, where Tg is found from the dilatometric tests, 
characterizes a measure of undercooling for a polymer. It is worth noting that "linear" 
Eqs. (31) are purely phenomenological and they cannot adequately describe the influence 
of temperature on the adjustable parameters in the close vicinities of critical points, Tg 
and Tc, where they result in corner points or even jumps of the graphs, whereas we suppose 
that the functions -B(T), C(T) and G{T) are rather smooth. 

5 Discussion 

Figures 5 and 6 demonstrate that the rate of increase, Fi, in the relaxation rate, F, 
with temperature reaches its maximum in the transition region (Tg,Tc), whereas in the 
sub-Tg and post-Tc domains this rate is essentially smaller. In other words, a decrease of 
temperature in the interval (Tg, Tc) by 1 K results in slowing down of the relaxation process 
which substantially exceeds that for temperatures out of this interval. This implies that 
changes in the internal structure of amorphous polymers between the glass transition 
point, Tg, and the critical point, Tc, are strongly affected by some additional micro- 
mechanism which does not "work" above the critical temperature and below the glass 
transition temperature. Bearing im mind that an amorphous polymer reveals a pattern 
of dynamic heterogeneity at evelated temperatures, whereas it exhibits a pattern of spatial 
inhomogeneity at low temperatures, the region (Tg,Tc) may be identified as a transition 
region from dynamic to static heterogeneity. Assuming this transition to be associated 
with the onset of surface energy of CRRs, we may conclude that an increase in the 
activation energy of amorphous polymers in the transition domain is driven by the growth 
of the surface energy of cooperatively relaxing regions. 

Figures 7 and 8 show that the rate of anelastic fiow, A, increases with temperature, 
T, in agreement with the theory of thermally activated processes. However, in contrast 
with Figures 5 and 6, these figures reveal different effects of temperature on the param- 
eter A: the rate of increase in A is maximal in the transition region (Tg,Tc) for PS, 



logB = Bo + BiAT, 



logC = Co + CiAT, 



logG = Go-GiAT, 



(31) 
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whereas the maximum of this quantity is reached in the post-Tc region for PMMA. Be- 
cause the Fisher temperatures, Tp, for these two polymers have not yet been determined 
experimentally, one can only speculate about an explanation for this observation. As 
two possible mechanisms for this phenomenon, we would mention (i) the effect of side 
groups in poly(methyl methacrylate) molecules, and (ii) the influence of parallelization of 
aromatic rings in polystyrene. 

Because the mobility of side groups in PMMA is rather high in the region (Tg, T^), they 
may create temporary junctions which strongly restrict mutual displacement of CRRs and 
slow down anelastic flow. The onset of temporary junctions is assumed to be driven by 
van der Waals forces and dipole interactions between side groups belonging to different 
chains. With the growth of temperature, T, the strength of these junctions decreases, 
whereas the intensity of thermal fluctuations increases, which implies that above the 
critical temperature, T^, the rate of anelastic flow in PMMA grows more rapidly, since 
the influence of temporary crosslinks between side groups becomes negligible (curves 2 
and 3 in Figure 8). 

The decrease in the slope of the graph A{T) for PS at temperatures above may 
be associated with parallelization of aromatic rings belonging to different macromolecules 
(inter-chain interactions) which implies (partial) molecular ordering and a decrease in the 
segmental mobility in the less cohesive space . Both these factors result in a decay in 
the rate of mutual displacement of CRRs. Below the critical temperature parallelization 
does not occur because the energy of thermal fluctuations is insufficient for an appropriate 
ordering of aromatic rings (curves 1 and 2 in Figure 7). 

Figures 9 and 10 demonstrate that the parameter B decreases below the critical tem- 
perature, Tc, reaches its minimum at T^, and afterwards, monotonically increases with 
temperature. It follows from Eq. (27) that B is inversely proportional to r, where r is 
the ratio of the rigidity of CRRs to the rigidity of LCDs. The increase in B (which is 
tantamount to a decrease in r) with temperature at T > Tc is quite natural, because 
the rigidity of less cohesive space grows (according to the theory of rubber elasticity, the 
elastic modulus is proportional to temperature), whereas the rigidity of CRRs decreases 
and vanishes at the Fisher temperature, Tp. The decrease in B (or the increase in r) 
with temperature, T, below means that under cooling of a sample, the rigidity of less 
cohesive space between CRRs grows faster than the rigidity of relaxing regions. This may 
be explained by the growth of the surface energy for CRRs, which substantially reduces 
mobility of strands that bridge rearranging domains and enhances their adsorption on 
the boundaries of CRRs, which, in turn, leads to an increase in the mechanical energy 
of chains in the less cohesive domains This explanation is supported by the fact 



that temperature-induced changes in the parameter B become substantially weaker in 
the sub-Tg region (Figure 10). 

According to Figures 11 and 12, the parameter C decreases with temperature in the 
sub-Tg region and increases in the post-Tg region. The decrease in C with temperature, T, 
below the glass transition point seems plausible, because it is equivalent to an increase in 
the rate of anelastic flow, see Eq. (24). The increase in C with temperature in the rubbery 
domain (T > Tg) is not rather pronounced (Figures 11 and 12 present semi-logarithmic 
plots of the data), and it may be explained by the onset of the dynamic heterogeneity 
in polymers: due to thermal fluctuations, CRRs disappear earlier than they perform 
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substantial displacements with respect to one another. This conclusion is confirmed by 
results depicted in Figure 11 which demonstrate that thermally induced changes in C 
practically vanish above the critical temperature, Tc, where dynamic heterogeneity plays 
the key role. 

Figures 13 and 14 show that the critical temperature, T^, cannot be identified from 
the dependence of the elastic modulus, G, on temperature: the parameter G monotoni- 
cally decreases with T in the entire post-Tg region. According to Figure 14, the rate of 
decrease in G is higher in the post-Tg domain than in the sub-Tg region. The latter seems 
physically plausible if we suppose that changes in G in the post-Tg region are affected 
by two processes: (i) thermally-driven decrease in the rigidity, c, and (ii) reduction of 
the surface energy for CRRs. A weak effect of the critical temperature on elastic moduli 
may be a reason why the effects studied in this work have not been previously revealed 
in conventional short-term tests. 

The temperatures, Tg and Tc, found as critical points for the dependences of material 
parameters on temperature, T, are collected in Table 2. The table demonstrates rather 
small scatters of data, which may serve as an indirect validation of the model. In par- 
ticular, it is found that the average glass transition temperature for PMMA found by 
matching observations in the mechanical tests is extremely close to that determined in 
the dilatometric tests (Tables 1 and 2). The latter conclusion is not surprising because 
the rate of cooling in dilatometric tests is comparable with the creep rate in long-term 
mechanical experiments. 

The critical temperatures for the two polymers are rather close to each other, which is 
in agreement with observations which evidence that PS and PMMA have similar values 
of most physical parameters [^. The latter conclusion is confirmed by the data listed in 
Table 3 which show quantitatively the same effects of temperature on adjustable param- 
eters for PS and PMMA (the constants which appear to be similar are printed in bold). 
The only exception from this rule is the rate of anelastic flow. A, which seems to depend 
strongly on the molecular weight and chemical structure of the polymers. 

6 Concluding remarks 

Experimental data are presented in long-term shear creep tests on polystyrene and poly- 
(methyl methacrylate) at various temperatures near the glass transition points. Constitu- 
tive equations are developed for the uniaxial mechanical response of amorphous polymers 
at small strains. The model is based on the hypothesis that the time-dependent behavior 
of polymers is governed by two micro-mechanisms: rearrangement of CRRs describes the 
viscoelasticity, whereas mutual displacement of CRRs reflects the anelastic phenomena. 
Adjustable parameters in the stress-strain relations are found as functions of tempera- 
ture. It is revealed that some critical temperature exists for the polymers, Tc ~ Tg -|- 20 
K, where the dependences of material constants on temperature substantially change. 
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List of captions 



Figure 1: The specific volume v cm^/g versus temperature T °C for polystyrene. 
Circles: observations. Solid lines: approximation of the experimental data by Eq. (1). 
Curve 1: vq = 0.9529, ao = 2.2213 x 10"^ K"^; curve 2: vq = 0.9173, ao = 6.4033 x 10"^ 

Figure 2: The specific volume v cm^/g versus temperature T °C for poly(methyl 
methacrylate). Circles: observations. Solid lines: approximation of the experimental 
data by Eq. (1). Curve 1: Vq = 0.8332, ao = 2.6172 x 10"^ K'^; curve 2: vq = 0.7965, 
tto = 7.3739 X 10-^ 

Figure 3: The creep compliance J GPa~^ versus time t s for polystyrene at a tem- 
perature T °C. Circles: experimental data. Solid lines: results of numerical simulation. 
Cmvc 1: T = 95.0; curve 2: T = 97.5; curve 3: T = 100.0; curve 4: T = 105.0; curve 
5: T = 110.0; curve 6: T = 115.0; curve 7: T = 120.0; curve 8: T = 130.0; curve 9: 
T = 140.0 

Figure 4: The creep complience J GPa~^ versus time t s for poly (methyl methacry- 
late) at a temperature T °C. Circles: experimental data. Sohd lines: results of numerical 
simulation. Curve 1: T = 92.5; curve 2: T = 97.5; curve 3: T = 100.0; curve 4: T = 102.5; 
curve 5: T = 105.0; curve 6: T = 107.5; curve 7: T = 110.0; curve 8: T = 112.0; curve 9: 
T = 115.0; curve 10: T = 120.0; curve 11: T = 125.0; curve 12: T = 130.0 

Figure 5: The relaxation rate F s^^ versus temperature T K for polystyrene. Circles: 
treatment of observations. Solid lines: approximation of the experimental data by Eq. 
(28). Curve 2: medium-temperature region, Fq = 50.73, Fi = 19373.5; curve 3: high- 
temperature region, Fq = 12.12, Fi = 4298.9 

Figure 6: The relaxation rate F s^^ versus temperature T K for poly(methyl metha- 
crylate). Circles: treatment of observations. Solid lines: approximation of the experimen- 
tal data by Eq. (28). Curve 1: low-temperature region, Fq = 21.35, Fi = 8759.6; curve 2: 
medium-temperature region, Fq = 60.10, Fi = 23245.7; curve 3: high-temperature region, 
Fo = 10.07, Fi = 3893.6 

Figure 7: The rate of anelastic flow A s^^ versus temperature T K for polystyrene. 
Circles: treatment of observations. Solid lines: approximation of the experimental data 
by Eq. (30). Curve 2: medium-temperature region, = 80.79, Ai — 30921.1; curve 3: 
high-temperature region, Aq — 37.52, Ai — 14371.1 

Figure 8: The rate of anelastic flow A versus temperature T K for poly(methyl 
methacrylate). Circles: treatment of observations. Solid lines: approximation of the 
experimental data by Eq. (30). Curve 1: low-temperature region, Aq = 23.11, Ai — 
9934.8; curve 2: medium-temperature region, Aq — 47.03, Ai — 18886.4; curve 3: high- 
temperature region, Aq — 95.73, Ai — 37851.8 

Figure 9: The dimensionless parameter B versus the degree of undercooling AT K 
for polystyrene. Circles: treatment of observations. Solid lines: approximation of the 
experimental data by Eq. (31). Curve 2: medium-temperature region, Bq = —0.4596, 
Bi = -0.0435; curve 3: high-temperature region, Bq = -3.3427, Bi = 0.0617 

Figure 10: The dimensionless parameter B versus the degree of undercooling AT 
K for poly (methyl methacrylate). Circles: treatment of observations. Solid lines: ap- 
proximation of the experimental data by Eq. (31). Curve 1: low-temperature region. 



18 



Bq = —0.3521, Bi = —0.0236; curve 2: medium-temperature region, Bq — —0.4778, 
Bi = -0.0904; curve 3: high-temperature region, Bq = -2.8531, Bi = 0.0789 

Figure 11: The parameter C GPa^ versus the degree of undercoohng AT K for 
polystyrene. Circles: treatment of observations. Solid lines: approximation of the ex- 
perimental data by Eq. (31). Curve 2: medium-temperature region, Cq — —6.6820, 
Ci = 0.1075; curve 3: high-temperature region, Cq = —4.4625, Ci = 0.0083 

Figure 12: The parameter C GPa^ versus the degree of undercooling AT K for 
poly (methyl methacrylate) . Circles: treatment of observations. Solid lines: approxi- 
mation of the experimental data by Eq. (31). Curve 1: glassy region, Co = —3.9026, 
Ci = -0.4004; curve 2: rubbery region, Co = -6.1357, Ci = 0.0908 

Figure 13: The elastic modulus C GPa versus the degree of undercooling AT K 
for polystyrene. Circles: treatment of observations. Solid line: approximation of the 
experimental data by Eq. (31) with Go = 0.1094 and Ci = 0.0699 

Figure 14: The elastic modulus C GPa versus the degree of undercooling AT K 
for poly (methyl methacrylate). Circles: treatment of observations. Solid lines: approx- 
imation of the experimental data by Eq. (31). Curve 1: glassy region, Cq = —0.2635, 
Ci = 0.0093; curve 2: rubbery region, Co = -0.1092, Ci = 0.0762 
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Table 1: Material constants determined by dilatometric measurements 



Polymer 


Tg °C 


"glass X 10^ 


"rubber X 10^ K'^ 


PS 


94.8 


2.12 


5.87 


PMMA 


99.5 


2.18 


5.87 



Table 2: The temperatures Tc and Tg determined as critical points for the dependences 
of adjustable parameters on temperature 



Parameter 


Tc °C (PS) 


Tc °C (PMMA) 


Tg °C (PMMA) 


r 


117.3 


113.6 


100.7 


A 


109.3 


116.3 


101.0 


B 


113.4 


112.9 


97.0 


C 


116.8 




103.5 


Go 






101.2 


Average 


114.2 


114.3 


100.7 



Table 3: Material constants determined by mechanical measurements 



Parameter 


Temperature 


PS 




PS 




PMMA 


PMMA 


r 


T<Te 




50.73 




19373.5 


To = 60.10 


Ti = 23245.7 




T>Te 


ro = 


12.12 


ri = 


4298.9 


To = 10.07 


Ti = 3893.6 


A 


T <T, 


Ao = 


80.79 


A,= 


30921.1 


Ao = 47.03 


Ai = 18886.4 




T>Te 


Ao = 


37.52 


A,= 


14371.1 


Aq = 95.73 


= 37851.8 


B 


T<Tc 


Bo = 


-0.46 


Bi = 


-0.04 


Bo = -0.48 


5i = -0.09 




T> 7; 


Bo = 


3.34 


Bi = 


0.06 


Bo = 2.85 


Bi = 0.08 


C 


T<Tc 


Co = 


-6.68 


Ci = 


0.11 


Co = -6.14 


Ci = 0.09 




T > Tc 


Co = 


-4.46 


Ci = 


0.01 


Co = -6.14 


Ci = 0.09 


G 


all 


Go = 


0.11 


Gi = 


0.07 


Go = -0.11 


Gi = 0.07 
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